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ABSTRACT:  In  this  paper,  the  matrix  pencil  method  has  been  utilized  for  estimating  the 

natural  resonances  from  different  transient  responses  recorded  along  multiple  look  directions  as 
a  function  of  time  after  the  incident  field  has  passed  the  structure.  The  novelty  of  this  article  is 
that  a  single  estimate  for  all  the  poles  are  done  utilizing  multiple  transient  waveforms  emanating 
from  the  structure  along  multiple  look  directions.  The  SEM  poles  are  independent  of  the  angle 
at  which  the  transient  response  is  recorded.  The  only  difference  between  the  various  waveforms 
are  that  the  residues  at  the  various  poles  are  of  different  magnitudes.  Some  of  the  residues  may 
even  be  zero  for  some  of  the  poles  indicating  that  the  contribution  from  certain  SEM  poles  may 
not  be  significant  along  that  look  direction.  Here  all  the  waveforms  are  utilized  providing  a 
single  estimate  for  the  poles  without  performing  an  arithmetic  mean  of  the  various  waveforms. 

1.  INTRODUCTION:  It  is  well  known  in  the  electromagnetics  literature  that  after  the 
incident  field  had  crossed  the  structure  of  interest,  the  time  domain  responses  can  be  modelled 
by  a  sum  of  complex  exponentials  [1-2].  In  the  Laplace  domain  this  is  equivalent  to  modelling 
the  transfer  function  of  the  system  by  the  poles  along  with  its  residues  or  in  terms  of  a  ratio  of 
two  polynomials  whose  roots  provide  the  poles  and  zeros  of  the  system. 

Many  methods  exist  in  the  published  literature  to  carry  out  such  a  parameterization  of  the 
source  free  transient  responses  of  the  system.  A  partial  survey  of  such  techniques  is  available 
in  [3,9,10],  Out  of  most  of  the  techniques,  the  matrix  pencil  method  had  proved  to  be  quite 
useful  [4,5]  because  of  its  low  sensitivity  to  background  noise  and  its  computational  ease  and 
efficiency. 

Now  when  the  transient  responses  from  the  object  of  interest  whose  SEM  poles  need  to 
be  found  out  is  looked  at  from  different  angles  both  in  azimuth  and  in  elevation,  the  residues  of 
the  poles  are  angle  dependent  whereas  the  SEM  poles  modeling  the  time  domain  waveforms  are 
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not  angle  dependent.  In  addition,  for  each  look  direction  there  are  two  possible  polarizations. 
One  could  also  use  both  polarizations  to  increase  the  number  of  waveforms  available. 
Conventionally,  to  estimate  the  SEM  poles  from  multiple  look  angle  data  one  takes  the  average 
of  all  the  various  look  directions  waveforms  and  then  obtains  a  single  waveform.  Then  a  sum 
of  complex  exponentials  is  used  to  fit  the  single  waveform  and  an  estimate  of  the  SEM  poles  is 
obtained  along  with  the  averaged  values  of  the  residues.  However,  this  is  not  a  good  approach 
if  the  signal-to-noise  ratios  of  the  different  waveforms  are  different  -  namely  in  some  of.  them 
the  transient  response  dies  down  quite  fast  whereas  in  some  of  the  responses  this  may  continue 
to  ring  for  a  long  time.  Hence  taking  an  average  of  those  two  classes  of  waveshapes  actually 
deteriorates  the  signal-to-noise  ratio  of  the  data.  This  is  because  by  taking  an  average  of  the 
signal  along  with  waveforms  where  the  signal  has  died  down  may  lead  to  an  unnecessary 
contamination  of  the  signal  by  noise.  In  this  paper  the  matrix  pencil  approach  is  applied  to 
obtain  a  single  estimate  for  the  SEM  poles  utilizing  simultaneously  all  the  transient  waveforms 
from  multiple  look  directions  and  without  averaging  them. 

In  the  next  section  the  matrix  pencil  method  is  presented  for  the  simultaneous  estimation 
of  all  the  SEM  poles  from  multiple  look  directions  without  averaging.  In  section  3,  the 
computational  procedure  utilizing  the  total  least  squares  singular-value  decomposition  based 
approach  is  presented  for  the  estimation  of  the  SEM  poles  from  multiple  look  directions.  This 
approach  has  been  found  to  be  most  robust  in  obtaining  estimates  for  the  poles  in  the  presence  ' 
of  random  noise  [6-8].  Section  4  provides  some  numerical  examples  utilizing  sample  simulated 
data  followed  by  conclusion  and  a  selected  set  of  references  where  additional  materials  are 
available. 

2.  Application  of  the  Matrix  Pencil  Method  for  Simultaneous  Estimation  of  the  SEM-Poles 

Utilizing  Waveforms  from  Multiple  Look  Directions 

Let  us  denote  the  transient  response  of  length  N+l  along  a  particular  look  direction  -k  by 
the  set  [YJ.  So  that  the  column  vector  [YJ  is  represented  by 
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[Ykv,)*,  -  ykO)5  ■  yk(N>]uo«.i)  (1) 

where  T  denotes  the  transpose  of  a  matrix.  The  elements  yk(j)  represents  the  values  of  the 
transient  response  at  the  jth  time  sample,  so  that 


yfc(j)  -  E  A^i)  exp[sj  AT]  ...  for  j  =0,  1,  2,  — ,  N  (2) 

i=t 

where  AT  is  the  sampling  time.  Each  transient  response  consists  of  the  same  M  SEM  poles  Sj 
which  are  to  be  solved  for  along  with  their  amplitudes  A^fi)  for  a  particular  look  direction  k. 
Please  note  that  M  is  also  an  unknown  along  with  the  SEM  poles  and  their  residues.  The  SEM 
poles  Sj  are  look  direction  independent  but  not  their  residues  Ak(i).  In  the  sampled  domain,  (2) 
can  be  rewritten  as 

yk(j)  =  E  A^i)  2jj  for  j  =0,  1,  2,  ....  N 

i-t 

where 

z.  =  exp[Sj  AT]  (4) 

It  has  been  further  assumed  that  all  the  waveforms  for  different  look  angles  k  =  l,  2, ...  K, 
have  been  sampled  uniformly  at  the  same  sampling  rate  AT  and  that  each  waveform  contains  the 


same  number  of  samples  N+l. 

'  Next  we  consider  two  matrices  [B,]  and  [BJ  defined  as; 
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[B,]s 


y,(0)  y2(0)  ...  yK(°) 

y,(i)  y20>  -  y*0) 

y,(N-l)  ya(N-l)  ...  yK(N-l) 


INxK 


(5) 


PM 


2JNxK 


y,0>  y20)  -  yKd) 

y,(2)  y2(2)  ...  yK(2) 
y,(N)  y2(N)  ...  yK(N) 


NxK 


(6) 


Now,  it  can  be  shown  that  the  two  matrices  [B,]  and  [B2]  can  be  decomposed  as  follows: 


~  PMnxM  PI  MxM  ^MxK 


(7) 


PMn*K  ”  JZJnxM  tZg]MKM  [A] 


MxK 


(8) 


where 


CZJ, 


MxM 


1  1 

h  *2 

Z*  *2 


N-l  N-I 


Z,  Zj 


M 

ZM 


Cl 


(9) 


INxM 


[I]MxM  =  is  a  diagonal  matrix  =  Identity  matrix 


(10) 


and  a  diagonal  matrix  [ZJ  containing  the  SEM  pole  set 
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[ZJ 


MxM 


z,  o 

h 

\ 

O  Z 


M 


MxM 


with 


tA] 


A,(l)  A,(l) 
A,(2)  Aj(2) 


|a,(M)  AjCM)  - 

Now  if  wc  consider  the  matrix  pencil 

[B,]  -  MB,] 


Ak(l) 

Ak(2) 

A^)jmxk 


(11) 


(12) 


(13) 


then  we  observe 

py  -  mb,]  -  [z,i  {izj  -  Mil}  [A]  (I4) 

% 

This  matrix  pencil  becomes  linearly  dependent  when  \  is  one  of  the  system  poles  as  then  the 
rank  of  ([ZJ  -  by  one  as  X  =  z.  Equation  (14)  can  be  transformed  into  a 

computationally  palatable  form  by  considering  the  ordinary  eigenvalue  problem  in  either  of  the 
following  forms: 

(15) 


PJ  [» ,r  -  m 


(16) 


m  -  X[B(]  [B J 

where  the  superscript  +  is  the  pseudo  inverse  of  the  respective  matrices.  The  pseudo  inverse  is 
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defined  in  terms  of  the  singular  value  decompositions  of  the  respective  matrix.  Let 


a  o 


4 


o 


[V,] 


H 


iJKxK 


(17) 


NxK 


=  [U,]  [£]  [V,]H 

where  (UJ  and  [V,]  are  two  orthogonal  matrixes,  i.e. 


[U,]-'  =  [U,]H  (18) 

[V,]-1  =  [V,]H  (19) 

where  the  superscript  H  denotes  the  conjugate  transpose  of  a  matrix.  Here  [X]  is  a  rectangular 
matrix  whose  diagonal  elements  are  related  to  the  singular  values  of  [B,].  In  summary,  we  have 
the  following  relationship 

(B, W  -  <*,  [(Uj]Nx,  for  c  =  1 , 2, ....  K  (20) 

and 


I^i^nxn  "  [Wnxi  :  Wnxi  :  —  :  ^un^]NxN  •••  (21) 

f^JxxK  -  [WkxI  :  l^Kx!  :  —  :  1Vk)rxi]KxK  "•  (22) 


Now  the  pseudo  inverse  of  [B,]  can  be  computed  from 
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[B.r  =  [V,]  [E]-1  [U,]H  (23) 

where 


(24) 

NxK 

It  is  interesting  to  observe  from  (7),  (8),  (15)  and  (16)  is  that  the  matrix  pencil  has  a  solution 
provided 

K  >  M  (25) 

i.e.  the  multiple  look  directions  must  be  greater  than  or  equal  to  the  number  of  poles  of  the 
system  to  be  estimated.  This  can  be  a  serious  limitation  in  many  cases  as  described  in  [11]  as 
the  number  of  SEM  poles  can  be  quite  large  for  many  practical  systems  and  it  may  not  be 
possible  to  provide  as  many  sensors  for  each  look  directions.  Hence  this  method  is  extended  to 
the  case  where  one  may  have  K  <  M.  If  K  <  M,  we  assume  that  N  »  K  or  M. 

To  deal  with  the  more  general  situation  we  consider  the  two  matrices  [D,]  and  [D2].  They 


[E]-1  = 


cl 


o 


o 


are  defined  by 
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[^J(L*l)xK-(N-L) 


y,(0)  y,(l)  -  y.(N-L-l)  y2(0)  y2(l)  -  y2(N-L-l) 

y,(l)  y,(2)  -  y,(N-L)  y2(l)  ya(2)  -  y2(N-L) 


y,(L)  y,(L+l)  ...  y,(N-l)  y2(L)  y2(L+l)  -  y2(N-l) 


and 


yK(°)  yK(0  -  y^N-L-1) 

yK0)  yK(2)  -yK(N-L) 


yK(L)  yK(L+i)  ...  yK(N-i) 


(27) 


(L*1)X[K-(N-L)1 


PJa, 


2J(Ul)xK-(N-L) 


y,(i)  y,(2)  -  y,(N-L)  y2(l)  y2(2)  -  y2(N-L) 

y,(2)  y,(3)  -  y.CN-L+1)  y2(2)  y2(3)  -  y2(N-L+l) 


y,(L+l)  y,(L+2)  ...  y,(N)  y2(L+l)  y2(L+2)  ...  y2(N) 


yK(l)  yK(2)  ...  yk(N-L) 

yK(2)  yK(3)  -  yK(N-L+l) 


yK(L+i)  yK(L+2) 


-  yK(N) 


i(L+l)x[K  *(N  -L)) 


(29) 


Next  it  can  be  shown  that  the  two  matrices  [D,]  and  [D2]  can  be  factored  into 


9 


P^H(L+l)xK-(N-L) 


=  [P] 


(L»l)xM  L  JMxM 


[UmxM  M 


Mx(L'M) 


[Q] 


(K-M)x[K-(N-L)] 


(30) 


and 


^2l(L*l)x[K(N-L)l  t^(L*l)xM  ^o-ImxM  [R^Mx(K  M)  [QJ(KM)x[K(N-L)] 


where 


[P]«, 


(L*l)xM 


11-1 


z,  z2 

2  2 
Z,  z2 


L  L 
Z,  Zz 


M 

2 

ZM 


(32) 


(L  +  l)xM 


[R] 


Mx(K-M) 


A,(l) 


O  =  A-(l) 


O  :  —  :  Ar(1)  O 

Aj(2)  :  A,(S)  :  -  :  AK(2) 

O  A,(M)  :  0  A^M)  :  :  O  AR(M) 


MxM 


(33) 


[Q] 


(K-M)x(K-(N-L)] 


[Q,] 


M  x(N=L) 


[Q,] 


1JM  x  (N-L) 


l(K-M)x[K-(N-L)] 


(34) 


and 
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[Q,] 


HMx(N-L) 


.  2  N-L-l 

1  Z,  Z,  -  Z, 

,  2  N-L-l 

1  z2  Z2  -  z2 


1  ZM  ZM 


N-L-l 


Mx(N-L) 


(35) 


In  addition  [I]  and  [ZJ  are  two  diagonal  matrices  which  have  been  defined  by  (10)  and  (11), 
respectively.  Now  if  we  consider  the  matrix  pencil 


[D2]  -  X[D.)  (36) 

and  their  equivalent  ordinary  eigenvalue  form  of  the  type 

[d2]  [D,r  -  m  w 

or 


[I]  -  X[D,)  [D 2r  (3g) 

then  when  A.  becomes  an  eigenvalue  of  either  (29)  or  (30),  and  its  value  is  equivalent  to  a  system 
pole. 

Once  all  the  poles  zu  i=l,  ....  M  have  been  computed  the  residues  at  the  poles  can  be 
computed  from  the  following  equation: 


y,(0) 

y2(0)  ■ 

••  yK(°) 

y,(D 

y20)  • 

••  yKd) 

y,(N) 

y2(N) 

yK(N) 

(39) 


l(N  +  l)xK 
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ri 

1  • 

••  1  ' 

r 

-i 

z,  z2  • 

ZM 

A,(l) 

A2(l)  - 

Ak(D 

2 

2 

2 

A,(2) 

A2(2)  - 

\(2) 

Z\  Zj  • 

■*  ZM 

: 

N  N 

7 .  7 ■ 

N 

•  •  7., 

A,(M) 

A2(M)  - 

Ak(M)_ 

(N  +  l)xM 

or  equivalently 


(40) 


[Y]  =  [Z]  •  [A] 


The  various  residues  can  now  be  computed  from  the  least  squares  solution  of  (31)  from 


[A]  =  [Z]*  [Y]  ={[Z]H[Z]}-'  [Z]H  [Y] 


(42) 


3.  COMPUTATION  OF  THE  SEM  POLES  UTILIZING  THE  TOTAL  LEAST  SQUARES 

In  order  to  deal  with  noisy  data,  the  formulation  of  the  previous  section  is  made  more 
robust  to  noise.  We  now  consider  the  composite  matrix  [D] 


y,(0) 

y,(D  •• 

•  y.(N-L-i) 

y.O) 

y,(2) 

-  y,(N-L) 

y,(L) 

y,(L+i) 

...  y,(N-l) 

y,(L+i) 

y,(L+2) 

-  y,(N) 

y2(°) 

y2(i)  - 

y2(N-L-l) 

y2d) 

y2(2)  - 

y2(N-L) 

y2(L) 

y2(L+i)  - 

y2(N-i) 

y2(L+i) 

y2(L+2)  - 

y2(N) 

[D]  = 
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yK(0)  yK0)  -  yK(N-L-i) 
yK(l)  yK(2)  -  yK(N-L) 


yK(L)  yK(L+i)  -  yK(N-i) 
yK(L+l)  yK(L+2)  -  yK(N) 


(43) 


l(L+2)x[K  -(N-L)] 


Please  note  that  [D,]  is  obtained  from  [D]  by  eliminating  the  last  row  and  [D2]  is  obtained  from 
[D]  by  eliminating  the  first  row.  We  now  perform  a  singular  value  decomposition  of  [D] 
according  to  (17)  as  follows: 


P"^(L*2)x(K-(N-L)l  [^](L-2)x(L*2)  t^(L+2)x [K -(N-L)]  ’  (N-l) x [K  (N-L)]  (44) 

To  combat  the  effects  of  noise  and  to  determine  the  order  M,  we  perform  a  singular  value 
filtering  of  [Z]  by  retaining  its  M  dominant  singular  values.  The  details  are  available  in  [4-8] 
and  is  omitted  in  this  paper.  Also  it  can  be  seen  that  the  ordinary  eigenvalue  problem  of 

[D2]  -  X[D,]  (45) 

can  be  transformed  into  the  following 


or 


t^2^(L+l)x  (L+2)  [^(L+2)x[K-(N-L)]  I^[K -(N-l)]  x  [K -(N-L)] 

^^“U(L+I)x(L+2)  [^](L»2)x[K'(N-L)]  t^][K -(N-l)]x[K -(N-L)] 


(46) 


[U2]  -  ^[U,] 


(47) 


where  [U2]  and  [U,]  of  (44)  are  obtained  from  [U]  by  eliminating  the  first  and  the  last  row, 


respectively. 


Then  the  poles  are  obtained  from  the  solution  of  either  one  of  the  following  four  ordinary 
eigenvalue  problem 

[U/  [U2]  -  MU/  [U,]  (48) 

[U2]  [U2]H  -  x  [U,]  [u/ 


[U2]  [U,]«  -  X[U,]  [U/  (49) 

[u,]H  [U2]  -  x  [U,]H  [U,] 

Once  the  poles  are  obtained  the  residues  at  the  poles  due  to  different  signals  measured  at  various 
look  directions  are  obtained  from  the  solution  of  (32)  through  the  use  of  (31). 


4.  NUMERICAL  RESULTS 

As  an  example  consider  a  square  plate  of  dimensions  lmxlm  (lying  in  the  x-y  plane) 
irradiated  by  an  electromagnetic  pulse  which  is  of  6 to  (light-meters)  in  duration  and  is  oriented 
along  the  E0  direction  with  magnitude  -377V/m.  We  are  observing  the  current  at  the  center  of 
the  plate.  The  waveforms  are  sampled  every  AT  =  0.1 1875  to  and  the  incident  pulse  dies  down 
after  12  to.  The  waveshape  for  the  y-directed  current  is  observed  after  130  time  samples  so  as 
to  ensure  that  the  incident  field  has  passed  the  metal  plate.  The  next  50  samples  are  taken  to 
estimate  the  SEM  poles  for  the  plate  by  observing  the  transient  field  arriving  from  different 
angles  of  incidence.  It  has  been  observed  that  only  5  poles  are  required  as  the  singular  values 
drop-off  beyond  10'8  in  evaluating  [D]  of  (43).  In  the  first  table  we  present  the  SEM  poles  along 
with  their  residues  for  seven  separate  incident  angles  of  0  and  (|). 
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Case  I;  0  =  0  and  <j>  =  90° 

S|  j,  =  -1.2  ±  jl.73 

A,  ,  =  ,45/±130° 

s3,4,  =  -.49  ±j.721 

A3  4  =  .37/118.5' 

s5  =  .703  A5  =  .0083 

In  this  case  the  incident  electric  field  Ee  of  amplitude  — 377V/m  is  impinging  on  the  plate  from 
0  =  0°  and  <(>  =  90°.  There  are  two  sets  of  complex  conjugate  poles  and  one  growing  exponential 
which  is  of  small  amplitude.  The  growing  exponential  is  non  physical  however,  it  is  the  error 


curve  fitting  of  the  data. 

Case  II;  0  =  10°  and  (j)  =  80° 

sli2>=  -1.35  ±jl.71 

A,  ,  =  .27/  +  38.2° 

s3  4  =  —.276  ±  j  1.241 

A,  ,  =  .19/+ 150° 

ss  =  .204 

> 

v> 

(1 

u> 

Case  III;  0  =  30°  and  <|)  =  50° 

s12j=  —1.06  ±  j  1 .75 

A,  2  =  .37/ +  32° 

s3  4  =  -.281  ±  jl.18 

A,  ^  =  .29/+ 143° 

s5  =  .041 

A5  =  .219 

Case  IV;  0  =  50°  and  <j>  =  75° 

s12  =  -.941  ±  j  1 .82 

A,  2  =  .44/ +  30° 

s3  4t  =  -.283  ±jl.l7 

A,  4  =  .36/ +138° 

ss  =  .07 

As  =  .325 
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Case  V;  0  =  20°  and  <j>  =  70° 


s,.2.  =  —1.57  ±  j  1 .49 

A,  ,  =  .9/T55.30 

s34>  =  -.278  ±  j  1.22 

A,  ,  =  .35/T1580 

s5  =  .14 

As  =  .454 

Case  VI;  0  =  10°  and  <|>  =  170° 

sl2>  =  —.963  ±  j4.ll 

A,  2  =  .145/=F55° 

s3  4t=  -.282  ±j  1.18 

A,  ,  =  .212/T 112° 

s5  =  .023 

A5  =  .208 

Case  VII;  0  =  30°  and  <)>  =  140° 

s)>2>=  —.999  ±  j2.42 

A,  2  =  .116/T310 

sM>  =  —.282  ±  j  1 . 19 

A,  4  =  .139/T  126° 

s5  =  .316 

As  =  .068 

As  it  can  be  seen  from  the  various  results  that  only  one  set  of  poles  around  —  .28+ jl.2  is  stable 
and  the  others  move  around.  The  various  poles  are  marked  in  Figure  1  through  the  various 
numericals  representing  the  seven  cases  described  above  by  roman  numerals. 

Next  we  utilize  all  the  seven  data  sets  to  obtain  a  single  estimate  for  the  poles.  Again  as 
before  5  poles  are  obtained.  This  single  estimate  of  the  poles  is  markes  by  **’  in  Figure  1. 
s1>2>  =  -1.13  ±  j  1.95 
s3  4_  =  -.28  ±  j  1.17 
s5  =  .039 


However  the  estimates  for  the  residues  are  different  for  different  look  angles. 


Case  I;  0  =  0°  and  (])  =  90° 
A1i2  =  ,35/t37.6° 

A3  4  =  ,273/t139° 

A5  =  .218 

Case  II;  0  =  10  and  <|)  =  80° 
A,  ,  =  .29/3=30.1° 

A3,4  =  .252/^131° 

A5  =  .21 1 

Case  HI;  0  =  30°  and  <|>  =  50° 
A,  2  =  .125/t45.2° 
A3>4  =  .152/t91.1° 

A5  =  .13 

Case  IV;  0  =  50°  and  <J)  =  70° 
A1i2  =  .126/t55° 

A3  4  =  .143/t86° 


A5  =  .117 
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Case  V;  6  =  20°  and  4  =  70° 
A1i2  =  .195/tI  1.7° 

A3  4  =  .213/^118° 

As  =  .186 

Case  VI;  0  =  10°  and  4>  =  170° 
A,  2  =  .077/t47.7° 

A,  4  =  .05/^152° 

A5  =  .036 


Case  VII;  0  =  30°  and  $  =  140° 

A,,  =  ,236/t46° 

A34  =  .159/^149° 

A5  =  .117 

From  the  theory  presented  in  this  paper  where  only  one  set  of  poles  has  been  estimated  using  all 
the  waveforms,  it  appears  that  there  is  a  possible  harmonic  relationship  between  the  two  decaying 
exponentials,  at  least  in  the  imaginary  part  and  as  expected  the  higher  frequency  is  damped  more. 
However  such  a  relationship  is  missing  when  the  poles  are  estimated  from  each  waveshape 


separately. 
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5.  CONCLUSION 

The  matrix  pencil  method  is  presented  for  estimating  the  SEM  poles  due  to  different 
transient  responses.  The  novelty  of  this  approach  is  that  only  one  set  of  the  SEM  poles  are 
estimated  from  the  multiple  waveforms.  However  the  residues  at  the  pole  sets  are  different  for 
different  waveshapes.  It  is  hoped  that  such  a  single  estimate  for  the  SEM  poles  will  be  more 
accurate  and  robust  to  different  orientation  of  polarizations  of  the  incident  fields  and  due  to 
various  effects  introduced  by  noise. 
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